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Abstract 

> 

Recent lattice simulations of (A$ 4 )4 theories in the broken phase show that : a) the 
shifted field propagator is well reproduced by the simple 2-parameter form 

at finite momenta but strongly differs for p — > b) the bare zero-momentum two- 

point function r 2 (0) = , f gives a value of = j^-r that increases when 

f ; VB=±VB 6 f r 2 (o) 

approaching the continuum limit. This supports theoretical expectations where vb is 
related by an infinite re-scaling to the 'physical Higgs condensate' Vr defined through 

d d v f = M\. New lattice data collected around the phase transition confirm 

this scenario. By denoting M SB = M h = 0(vr) the scale of the broken phase, 
our results suggest the existence of a 'hierarchy' of scales r 2 (0) <C M| B <C v% that 
become infinitely far in the continuum limit. This may open unexpected possibilities 
to reconcile an infinitesimal slope of the effective potential with finite values of Mh 
and accomodate very different mass scales in the framework of a spontaneously 
broken theory. 



1 Introduction 



Spontaneous symmetry breaking through an elementary scalar field is the basic 
ingredient for the origin of particle masses in the Standard Model of electroweak 
interactions. Traditionally, the 'condensation' of a scalar field, i.e. the transition 
from a symmetric phase where (<&) = to the physical vacuum where ($) 7^ 0, 
has been described as an essentially classical phenomenon in terms of a classical 
potential ('B=Bare', A > 0) 

with non-trivial absolute minima for a constant value = ±vb 7^ 0. In this 
picture, by expanding around the absolute minima of the classical potential say 
$_b(x) = vb + h(x), one predicts a simple relation 

Ml = V£(*b = vb) (2) 
between the 'Higgs mass' and the quadratic shape of the potential at the minima. 

In the quantum theory, and on the basis of perturbation theory, Eq.fl2|) is be- 
lieved to represent a good approximation by simply replacing the classical potential 
with the quantum effective potential V e s- However, beyond perturbation theory, 
there is an alternative description of spontaneous symmetry breaking 0-44I where 
M% and the curvature of the effective potential at the non-trivial minima are dif- 
ferent physical quantities related by an infinite renormalization in the continuum 
limit of quantum field theory, with potentially important consequences for parti- 
cle physics and cosmology. The above result has been obtained from gaussian and 
post-gaussian|J approximations to the effective potential. In view of the relevance 
of the issue and for the convenience of the reader, we shall briefly recapitulate the 
main result in the simplest case of the gaussian approximation to the energy density 
^gVpbi^A !]|" 10 1- This is obtained from the expectation value of the hamiltonian 



in the class of trial states with a constant = <Pb and shifted-field euclidean 

propagator 

with a variational mass parameter Q. By minimization one gets the coupled equa- 
tions 

n 2 = r + ^ + ±i (n) (4) 



and 

> 

= r„LOu A — 

6 2 



= r 0(pB + ^ + ^1,(0) (5) 



1 



with 



d 4 p 



(2vr) 4 p 2 + n 2 



(6) 



By using Eq.(f|), one can define a <£>b — dependent mass Q = Q((ps) and the gaus- 
sian approximation to the effective potential Vg((Pb) = £g[<Pb, ^(^b)] whose first 
derivative has the simple form M 



dV< 



G .. rrv2, . , Vbi 



(p B [tt (tp B ) - 



(7) 



In this way, spontaneous symmetry breaking is associated with those absolute min- 
ima ifs = ±^b 7^ where 



n 2 (v B ) 



Now, the zero-momentum two-point function (the inverse susceptibility) defines the 
quadratic shape of the potential at the non-trivial minima 



X 



T 2 (P = 0) 



d 2 V r 



cff 



dip 2 B 



Ml 



Vb=±v b 



where we have introduced a re-scaling factor Z = Z v defined through Eq. ( 
In the gaussian approximation, by using the identities || 

j o (0)-/ o (n) = y[/_i(n) 1 1 



7T- 2 



and the definition 



we obtain 



dlo 
dQ 



d 2 V ( 



G 



dip 



n 2 (<p B ) - 



B 



\<P 2 B _A/-i(«) 
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(10) 



(12) 



In the symmetric vacuum where ip B = 0, one finds 



<PV, 



G 



d(p 2 B 



n 2 (o) 



ip B =0 



while, in the broken-symmetry phase, the alternative expression 



d 2 V G 
dip% 



M, 



y B =±v B 



-f/-i(M fe ) 
h l + \l_ x {M h ) 



(13) 



(14) 
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For an evaluation of Eq. (|T4|) , it is convenient to re-define the bare mass term as 



r o = -^ o (0)+A (15) 



so that the gap-equation at the absolute minima reduces to 



= 4 + |-^[L l( J« + l| (16) 



( A = corresponds to the 'Coleman- Weinberg regime ' (12| where f2(0) = 0). By 
defining 

'-W 3 !?^ (17) 

and /i 2 = A 2 e _i ^, A = ^y, Mf t = oj 2 /i 2 , Eq.(|l|) becomes 

u 2 \n^ = y (18) 

Finally, by replacing into Eq.([b|), we obtain the gaussian approximation prediction 
for Z^ in the broken phase 

J_ = r 2 ( P = o) = i^in^ 2 

Z v Ml |-^ln^ 

that (in the range > uj 2 > 1) amounts to 



24tt 2 1 , . 

Eq.([20|) can be considered from two qualitatively different points of view. On one 
hand, for given values of the dimensionless parameters y and uj 2 , it implies a be- 
haviour Z v ~ \. Thus, in a continuum limit where A — > oo and fi 2 (and hence 
M 2 = /i 2 c<j 2 ) is kept fixed, one gets Z v ~ In — > oo implying that the 'Higgs mass' 

ft. 

M/j and the curvature of the effective potential are different physical quantities. In- 
deed, they would be related by an infinite renormalization in the continuum limit 
of quantum field theory. 

On the other hand, for any given A, i.e. when /i 2 and A are kept in a fixed 
ratio, Eq. (|20|) predicts that Z v should also become larger while decreasing uj 2 . This 
corresponds to increase y from negative values trying to approach the minimum value 
uj 2 = 1. However, in this latter case, Z v will not diverge. In fact, before reaching the 
value u 2 = 1, the gaussian effective potential exhibits a (weakly) first-order phase 
transition to the symmetric phase M. 
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We stress that the gaussian-approximation prediction of a non trivial has 
been tested with precise lattice simulations [T3|, TA\ performed in the Ising limit of 



the theory. The data show that by approaching the critical value of the hopping 
parameter k c ~ 0.0748 |15] from the broken phase, the quantity = M\x rapidly 



increases above unity. Thus, in the broken phase, Z^ is different from the more 
conventional quantity Z = Z prop associated with the two-parameter form for the 
shifted-field propagator 

= (21) 



According to 'triviality' [16|, this latter quantity should exhibit a continuum limit 
-Zprop — > 1 for consistency with the Kallen-Lehmann decomposition. The point is 
that, in the broken phase, the data deviate from Eq.(]2"l"D when p — > presaging an 
even more dramatic difference at p = 0. Furthermore, the discrepancy between Z v 
and Z prop becomes larger in the limit k — > k c . Therefore, it cannot be explained 
by introducing residual perturbative corrections that, according to 'triviality 1 JT^] . 



should become smaller and smaller in the same limit. On the other hand, in the 
symmetric phase where Eq.(|2T|) describes the lattice data down to p = |T3J], one 
gets Z v = Z prop ~ 1 as expected on the basis of Eq . (|i~3D . 

After this general Introduction, we shall report in this Letter further numerical 
evidence for the validity of this theoretical framework by studying the uj— dependence 
of Eq. fl20|) . Strictly speaking, to this end, one should work at fixed A by changing 
the bare mass in the full two-parameter A$ 4 theory, the Ising limit being just a one- 
parameter model. However, even in the Ising limit we can reproduce the situation 
of a two-parameter theory if we perform a finite-temperature simulation. Namely, 
by changing the value of k in the broken phase, we know from refs.fll3|,|i~4| that the 
zero-temperature value of Z v increases when k, — > k c . However, for any fixed value 
k > n c , by increasing the temperature, we shall approach the phase transition. In 
this case, we expect the lattice data to exhibit a qualitatively different increase of 
Z v , analogously to the effect induced in Eq.(|^) by approaching the phase transition 
at fixed A in the full two-parameter theory. 



2 Finite-temperature lattice simulation 



The transition to finite-temperature field theory is easily performed by following the 
formalism introduced in refs.|l7|]. In the case of bosonic fields the partition function 
is 



Z{P) 



£>$ e -Io dT I d3xC E 



(22) 
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where (3 = 1/T is the inverse temperature, He is the Euclidean lagrangian density 
and the functional integral is performed over all periodic field configurations with 
$(x, r = 0) = $(x, r = (3). In Fourier space, the compactification of the time 
direction leads to a discrete energy spectrum with energies u) n = ^hl and, therefore, 
to the form of a free propagator 



1 



p 2 + pV + M 2 



(23) 



As it is well known |]18| , a Monte-Carlo simulation with periodic boundary con- 
ditions on an asymmetric x L t lattice is equivalent to a finite-temperature calcu- 
lation at (3 = L t = 1/T. To this end, a one-component (A$ 4 )4 theory 



\ + e,) - $(x)) 2 + T ^\x) + ^ 



x 



is conveniently studied in the Ising limit 

Rising = _K + ^)^( X ) + <K X _ e M )0(x)] 



(24) 



(25) 



with $(x) = v^2k0(x) and where 0(x) takes only the values +1 or —1. 

To perform Monte-Carlo simulations of this Ising action, we have used the 



Swendsen-Wang |19| cluster algorithm. Statistical errors can be estimated through a 
direct evaluation of the integrated autocorrelation time [BO], or by using the "block- 



ing" jn|] or the "grouped jackknife" [22] algorithms. We have checked that applying 
these three different methods we get consistent results. Lattice observables include: 

(i) the bare magnetization, vb = (\&\), where <3> = ^2 x ^(x)/ L 4 is the average 
field for each lattice configuration. The broken phase is found for k > k c where 
k c ~ 0.0748 |i~5| in an infinite volume. For any k > k c , we expect to detect a phase 
transition to the symmetric phase at some sufficiently small value of L t . 



(ii) the zero-momentum susceptibility 

X = L* [<|$| 2 )-(|$D 2 ] 



(26) 



(iii) the shifted-field propagator 



G{p) = I ^exp(ipx)($(x) -v B )($(0) -v B ) 



(27) 



where p = (p4,p) with p = and p 4 = j^n^ where (^,112, 11,3,714) are integer- 
valued components, not all zero. 
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To extract the 'Higgs mass' one has to preliminarly compare the lattice data 
for G(p) with the 2-parameter formula 

Gg(n,p) = - ■ Z T P 2 (28) 



p + P4 + m 



latt 



where m\ att is the mass in lattice units and p fJL = 2 sin 

Obviously, by setting L t = L s we re-obtain the zero temperature case studied in 



refs.[]13|)p[. We now repeat the basic point of |L4|. To realize how good the fit with 



Eq.fl28|) can be, we start at T = in the symmetric phase at k = 0.0740 on a 20 4 



lattice. In Fig.l, we report the same data of ref . | Hfl where the scalar propagator 
has been suitably re-scaled to show the very good quality of the fit to Eq. (p8"|) . The 
value of Z prop is indicated by the dashed line while Z v is reported in Fig.l as a black 
point at p — 0. Notice the perfect agreement between Z v and Z prop . 

We now select a value k = 0.07512, in the broken phase, where the 2-parameter 
fit to the propagator data yields a comparable value for m latt by using a 32 4 lattice. 
However, unlike Fig. 1, the fit to Eq. (^8l), though excellent at higher momenta, does 
not reproduce the lattice data down to zero-momentum. Therefore, in the broken 
phase, a meaningful determination of Z prop and mi att requires excluding the lowest 
momentum points from the fit. The lattice data are shown in Fig.2. 

By fixing k, = 0.07512, we have now performed computations on asymmetric L 3 x 
L t lattices to simulate a finite-temperature system. Notice that well past the phase 
transition, i.e. for very small L t , consistency requires ap — - > limit of the propagator 
such that Zy and Z prop agree to good accuracy. In fact, if running on smaller lattices 
one is able to re-obtain the same conditions of the symmetric-phase calculation in 
Fig.l, we can exclude the presence of finite-size artifacts in the zero-temperature 
simulation and interpret the discrepancy between Z^ and Z pTop as a real physical 
effect due to the presence of the scalar condensate. At the same time, the phase 
transition associated with spontaneous symmetry breaking can be understood as a 
real condensation process where the Z v — Z prop discrepancy represents a distinctive 
non-perturbative feature. 

This expectation can be checked in Fig. 3 by comparing the data taken on the 
32 3 x 2 lattice with those taken on the 32-times bigger 64 3 x 8 lattice when the 
system, however, is still in the broken phase and a meaningful determination of 
mi a tt through Eq. (|28|) requires excluding the lowest momentum points from the fit. 



The data for the physical observables from different lattice sizes are reported in 
Table 1. For all lattice sizes there is a clear evidence for a phase transition in the 
region 6 < L t < 8 where the system crosses from the broken into the symmetric 
phase (see Fig. 4). 

We stress that the conventional interpretation of 'triviality', assuming Z v = 
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Z pr0 p ~ 1, predicts that approaching the phase transition one finds m 2 — > and 
X — > oo in such a way that m 2 x remains constant. However, our data for the 
quantity Z v = m 2 att x, reported in Fig. 5, exhibit a clear increase in the critical 
region while the corresponding value for Z prop remains remarkably constant (see 
Table 1). This confirms once more that Z v , as defined in Eq.(|9]), and Z piop , as 
defined in Eq . (|2l|) , are different physical quantities. 



3 Summary and outlook 

The finite-temperature simulations presented in this paper confirm the results of 
refs.[rT3|,|n|. Both provide strong numerical evidence that, close to the continuum 
limit of quantum field theory, a naive perturbative description of a spontaneously 
broken (A$ 4 )4 is missing very important effects. Although the theory is 'trivial' 
there are non-perturbative collective effects at (and near) zero momentum produc- 
ing the observed large deviations from the standard free-field-like form of the prop- 
agator Eq. (|28|) . These deviations are responsible for the non-trivial re-scaling factor 
Z^ = Mix that increases rapidly when k — > k c in the zero-temperature case and 
when approaching the phase transition in the finite-temperature simulations. This 
provides additional numerical evidence for the validity of Eq.(^0]) that represents a 
distinctive prediction and can be used, for instance, to reconcile a finite value of 
Mh with an infinitesimal quadratic shape of the effective potential at its minima. 
Indeed, with a divergent Z v , the usual definition of the physical vacuum field vr 



d 2 V c 



elf 



= M 2 (29) 

<p R =±v R 



is related by a non-trivial infinite re-scaling to the bare 'Higgs condensate' vb in 
Eq.([T^). Therefore, although vb diverges in units of Mh, one can obtain a continuum 
limit where both and vr are finite quantities setting the scale of the spontaneously 
broken phase Msb = = 0(vr). In this sense, a divergent Z^ means that 
spontaneous symmetry breaking introduces a mass 'hierarchy' [|J] where T 2 {p = 
0) < M| B < v 2 B are infinitely far scales in the continuum limit since 

M 2 i) 2 A 

^ = F7 ^T = ^~ln^-oo (30) 
T 2 (p = 0) v R M h 

Another interesting point concerns the actual form of the energy spectrum for the 
long- wavelength excitations of the broken phase. In fact, the observed differences of 
the propagator from Eq. (f2T|) for p — > imply that the energy spectrum of the broken 



phase E(p) sizeably differs from ^p 2 + M| when p — >• 0. In this respect, the results 
of ref.[14j] show that the spectrum E(p) approaches the form a/p 2 + M\ only at large 
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p 2 . Also, E(0) < M/j and their difference increases in the continuum limit. Such a 
difference, detected by studying the time-slices of the connected correlator at various 
values of p, has no counterpart in the symmetric phase. In this latter case, the form 
E(p) = a/p 2 + m 2 is found |14| to reproduce the lattice data to high accuracy for all 
p down to p = so that, in this case, one gets -E'(O) = m. A particularly important 
question concerns the stability of the results for the energy-gap E(0) in the broken 
phase and for the 'Higgs mass' Mh controlling the higher-momentum behaviour of 
the propagator. The results of |l4j for k = 0.076, the case studied by Jansen et al. 
[ p4| on a 20 4 lattice, show that the value of Mh, extracted from the set of higher- 
momentum data where one gets a good fit with Eq. (pT|) , is remarkably stable for 
variations of the lattice size from 20 4 to 32 4 . New preliminary data [^5| show that 
the energy spectrum E(p) remains also stable, at least for not too small values of 
p 2 . On the other hand, by increasing the lattice size up to 32 4 , our present data 
suggest a decrease of E(0) that, therefore, differently from Mh, may represent an 
infrared-sensitive quantity. A complete discussion of this point requires, however, 
more statistics and will be presented in a forthcoming paper [p5l . 
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lattice #configs. 1/L t < |$| > m latt x ^ z P ro P 

32 3 x 2 700K 0.5 0.015980(42) 0.494362(917) 26.25(0.14) 0.963675(6389) 0.960179(1167) 

32 3 x 3 835K 0.333333 0.024602(47) 0.261363(681) 93.28(0.34) 0.957328(6074) 0.955672(898) 

32 3 x 4 650K 0.25 0.035326(81)) 0.180635(6834) 253.60(1.09) 1.243208(94220) 0.959723(1649) 

32 3 x 5 875K 0.2 0.050104(105) 0.139061(5478) 614.61(2.36) 1.785661(140851) 0.959713(1306) 

32 3 x 6 725K 0.166667 0.071158(285) 0.125006(6219) 1376.14(9.58) 3.230797(322248) 0.959374(1373) 

32 3 x 8 725K 0.125 0.117674(375) 0.171012(7257) 530.55(3.85) 2.331141(198569) 0.959329(1677) 

32 3 x 32 400K 0.161778(130) 0.20623(409) 193.1(1.7) 1.233876(501) 0.9551(21) 

48 3 x 3 700K 0.333333 0.013413(27) 0.261399(594) 93.65(0.35) 0.961384(5645) 0.955420(834) 

48 3 x 4 460K 0.25 0.019187(47) 0.165589(2391) 255.35(1.19) 1.051909(30770) 0.956229(1257) 

48 3 x 5 635K 0.2 0.027706(72) 0.118999(2923) 655.07(3.26) 1.393679(68817) 0.958973(1088) 

48 3 x 6 400K 0.166667 0.044522(193) 0.114195(6459) 1919.71(15.16) 3.761108(426501) 0.960837(1577) 

48 3 x 8 170K 0.125 0.116546(516) 0.157074(7182) 826.14(20.94) 3.062311(290602) 0.958064(2392) 

48 3 x 9 185K 0.111111 0.136291(415) 0.180302(6762) 464.83(11.35) 2.270273(179082) 0.955628(2456) 

64 3 x 2 390K 0.5 0.005626(18) 0.494478(1463) 26.14(0.17) 0.9604047(8320) 0.959913(3017) 

64 3 x 3 250K 0.333333 0.008712(31) 0.260004(681) 93.52(0.74) 0.9498651(9044) 0.952527(1263) 

64 3 x 4 480K 0.25 0.012439(33) 0.163052(2724) 259.16(4.95) 1.035163(39846) 0.954254(1341) 

64 3 x 5 200K 0.2 0.017980(89) 0.110235(4120) 676.43(19.80) 1.234948(99135) 0.955591(1730) 

64 3 x 6 90K 0.166667 0.029673(263) 0.090056(9592) 2105.82(24.60) 2.565853(547407) 0.957225(2888) 

64 3 x 8 240K 0.125 0.118429(299) 0.162401(7698) 803.64(17.45) 3.184376(309700) 0.957547(2466) 

64 3 x 9 245K 0.111111 0.136991(178) 0.177897(4993) 417.49(7.65) 1.985071(117215) 0.955239(2092) 

64 3 x 10 150K 0.1 0.1462464(180) 0.180102(4437) 311.83(4.56) 1.519625(78097) 0.952393(2373) 

64 3 x 14 80K 0.0714286 0.158104(137) 0.193444(5602) 221.91(3.73) 1.247618(75245) 0.948202(3055) 



Table 1: Summary of the simulation runs for the 4d \tp A model at finite temperature (Ising limit, k = 0.07512). 
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Figure 1: The lattice data for the re-scaled propagator at k — 0.0740 in the sym- 
metric phase. The zero-momentum full point is defined as Z v = m^ tt x- The dashed 
line indicates the value of Z r _ 
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FIGURE 2 
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Figure 2: The lattice data for the re-scaled propagator at n — 0.07512. The zero- 
momentum full point is defined as = Tn^ tt x- The very low momentum region is 
shown in the inset. The dashed line indicates the value of Z DTOD . 
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FIGURE 3 
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Figure 3: The same as in Fig. 2 for L s = 32 and L t = 2, and L s = 64 and L t = 8. The 
two sets of data are normalized in terms of the corresponding lattice masses shown in 
Table 1. The dashed lines indicate the same zero-temperature value Z prop = 0.9551 
as in Fig.2 
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FIGURE 5 
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Figure 5: The lattice data for Z v = mf att x for L s = 32,48,64 vs. the temperature. 
The dashed line indicates the zero-temperature value Z prop = 0.9551 as in Fig.2 
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